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We examine the feasibility of enhancing the fundamental radiative interactions 
between distant atoms. We present general arguments for producing 
enhancement. In particular, we show how giant dipole-dipole interaction can be 
produced by considering dipoles placed close to micron sized silica spheres. The 
giant interaction arises as the whispering gallery modes can resonantly couple 
the dipoles. 

Introduction 

It is well known from the early days of quantum electrodynamics^ that the 
exchange of a photon between two atoms leads to the radiative coupling between 
the atoms. It is also known that the strength of such a radiative coupling is 
significant only if the distance between the two atoms is much smaller than the 
wavelength of the transition involved in the coupling.. In this work we 
demonstrate the possibihty of remarkably large radiative coupling between atoms 
which could be separated by distances much larger than the wavelength. In what 
follows, we concentrate on one specific type of the radiative coupling, viz., the 
retarded dipole-dipole interaction. We begin by recalling how the radiative 
coupHng arises. 

Consider, for example, two atoms A and B located at positions and , 
respectively. Let at time ? = the atom A be in the excited state , and the 
atom B in the ground state | ) . Such a selective excitation can be achieved by 
using an optical fiber tip. We assume an allowed dipole transition between 
and I ^a) • Clearly we can have a transfer of excitation from the atom A to B, i.e., 
the initial state of the system has changed from Ica'^b) to I^A'^b)- The 
electromagnetic field is in vacuum state at ? = and is again in vacuum state 
after the transfer of excitation is completed. This process \ej^,gg^ to |^A'^b) 
mediated by the vacuum of the electromagnetic field. This transfer process can 
also be described by an effective interaction, which depends only on the atomic 
degrees of freedom. Thus, the atoms A and B get coupled via the radiative 
interaction. The amplitude for such a process can be calculated in the framework 
of quantum electrodynamics using the second order perturbation theory. Let d„ 

be the dipole matrix element for the transition | e„ ) ^ | ^„ ) • As shown in the 
appendix a quantum electrodynamic calculation shows that the effective 



interaction is related to the dyadic Green's function Gij{r^,rg,(o) associated with 

the wave equation for the electromagnetic field. The dyadic Green's function is 
related to the field (r^ , CO) produced by an oscillating dipole at the position 

E,{?^,CD) = ^Y,Gy{?A,?j„CO)pj . (1) 



The function G is well known to be' 
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The role of G in determining the radiative interactions can be easily appreciated 
by examining the classical equation of motion for a dipole at the position 
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Here E{rj^,t) is the net field at location A and it should also include the self- 
fields^, i.e., its Fourier component is given by 

E{rA,Co)^^ [G(r^ , , CO).Pa + Gir^ , rg , CO).Ps ] . 



(4) 



We will now make slowly varying envelope approximation 
p,(t) = pJ,e-''^'+c.c., (5) 
where H„ represents the orientation of the dipole at . We can then obtain the 
reduced equation for p 
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Thus the net complex coupUng between two dipoles is given by 
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The real and imaginary parts of the Green's tensor determine the strength of the 
radiative interaction. The radiative interaction is an oscillatory function of 

i^A~^B) is especially important if — - « 2;r . If the distance 

c 

between the two dipoles is much larger than a wavelength, then the radiative 
coupling is negligibly small. The question that we ask is - is it possible to 
enhance the radiative coupling between atoms which are separated by a distance 
that is many times the wavelength of the atomic transition. As noted above the 
radiative interaction arises from the interaction with the vacuum of the 
electromagnetic field. Thus we should be able to enhance v by suitably 

modifying the vacuum. In other words, if the structure of G for the modified 

vacuum were such that G had a resonance at frequency corresponding to the 



atomic transition frequency, then v can be enhanced by a large factor. The actual 
value of the enhancement factor will depend on the widths of the resonances in 

G . There have been several studies on surface enhanced dipole-dipole 
interaction A situation which is particularly attractive is shown in Fig.l., 
where we place two atoms near the surface of a micro-sphere with dielectric 
constant e . It is known that micro-spheres support whispering gallery modes 
(WGM) and some of these modes with frequency say (Oq can have ultra-narrow 

width. In other words, if ftJ^ is in the vicinity of (Oq , then G has primarily a 

resonant character and a small background. In what follows, we demonstrate a 
very large enhancement of the radiative interaction between two atoms separated 
by a distance of, say, two orders of magnitude larger than the wavelength of 
transition. Note that in recent times whispering gallery modes ^ (WGM) have 
attracted considerable attention following the work of Braginsky et al who 
showed the possibihty of WGM's with quahty factors Q>10^, effective 
quantization volume of the order of 10'^ cm"^. These two properties make WGM's 
specially attractive in many different applications in nonlinear optics and in 
quantum electrodynamics Braginsky et al had already reported the possibihty 
of extremely low threshold for optical bistability. Haroche and coworkers have 
extensively studied the properties of WGM's These authors even developed 
newer techniques for studying such modes. Arnold and coworkers studied 
WGM induced enhancement of energy transfer, albeit for smaller spheres (of size 
= 5iMn ), though as will be seen in our work that much larger enhancement can 

be achieved by using larger spheres. Kimble et al and the present authors 
considered interaction of WGM's with atoms. In particular, Dutta Gupta and 
Agarwal considered micron sized spheres doped uniformly with resonant 
atoms. They demonstrated how effects hke vacuum field Rabi sphttings can be 
studied using WGM's. The lasing characteristics of doped micro-spheres have 
also been studied extensively 

Structure of G 

In order to obtain G for the situation shown in Fig.l, we have to solve Maxwell's 
equations subject to the boundary conditions at the surface of the sphere. The 

Greens function consists of two contributions " - (a) G^^^ - which is the same as 

in free space, i.e., in the absence of the micro-sphere; (b) G^'^ - which depends 
on the presence of the micro-sphere. The basic equations for these Green's 
functions are 

2 

V X V X G*"^ (F, r, (o)-^ eM^°^ {r, r', a) = 47a5{? - r), (8) 

c 

Vx'VxG^'^ir,r,co)-^eoG^'^ir,r,co)=0 . (9) 
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These Green's functions are given in terms of the vector spherical harmonics as 



G*"^ (r,, rs,co) = ik^ > > (2 - 5„,o ) \ , , , , , 
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In Eqs.(lO) and (11) (r,) represent the greater (smaller) of and r^, 
kQ=(ol c . The vector spherical harmonics M and iV are given by 
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The vector functions M^^^ and A'^*^^ are obtained from Eq.( 12) by replacing 
by /j^^^ . A and B in Eq.(ll) are determined by the following expressions 
^ in (ka)[ko^jn ikoa)Y-j„ (kQa)[kaj„ (ka)y 



ftf ikoa)[kaj„ (ka)Y-j„ {ka)[k^ah';^^ {k^a)\ 
£j„ (ka)[k^^aj„ (k^^a)]'-j„ {k^)a)[kai„ {ka)\' 



(13) 
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/z*'' {k^^a)[kaj„ (ka)]'-ej„ (ka)[k^^ahl^> (k^^a)]' 

In (13) and (14) k^k^^e . The vanishing of the denominators in (13) and (14) 
yields the various modes of the system 

ft« {k^a)[kai, (fca)]'-7„ {ka){k^ahf^ {k^a)\= , (15) 

ft«(^o«)[%'„(^)]'-e/;(^«)[^o«/^,*"(^o«)]'= . (16) 
Depending on whether Eq(15) or (16) is satisfied, the modes are labeled as a- 
(TM) or b- (TE), respectively, since the radial component of the magnetic or 
electric field is absent in the corresponding field distribution. The roots of 
Eqs.(15) or (16) and hence the modes can be labeled by two integers n and /, 
where n gives the mode number and / gives the order number. Physically / 



determines tiie number of extrema in the radial distribution of the field. The 
smaller the value of I, the more the field is concentrated near the inner edge of 
the sphere. In the ray picture this corresponds to a ray, confined near the inner 
edge, undergoing multiple total internal reflections from inside the micro-sphere, 
justifying the name "whispering gallery mode". Under plane wave illumination 
(as in extinction studies) a variely of such modes can be excited. Especially 
interesting are the modes of larger spheres with large n, for which the 
corresponding quahty factors can achieve dramatic values (order of 10^"^" in the 
visible domain). Note that it is extremely difficult to achieve such figures with 
conventional Fabry-Perot cavities. 

It is now clear that at chosen discrete set of mode frequencies satisfying Eqs.(I5), 

(16), the poles in A or B can lead to dramatic enhancements in G**^ and thus in 
radiative interactions. In order to demonstrate the enhancement and to facilitate a 
comparison of small and large spheres to reveal the potentials of the larger 
system, we picked sihca spheres of two radii, namely, 200 |im and 5 |im. For the 

larger (smaller) sphere out of the variety of modes we selected a23i2,i67 (^39,1) 
mode with frequency (0^= 1.8026933000486 |xm' (1.096835618 |xm') and 

quality factor Q = 1.28x10^ (1.22x10^). In the next Section we define the 
quantities of interest and present numerical results pertaining to these modes. 

Giant Enhancement in the Radiative Interaction 

We now discuss the relevant quantities capturing the giant enhancement in 
radiative interactions. We also note the presence of the self-terms in the equation 
of motion for the dipoles. Such self-terms give rise to the usual lifetimes and 
frequency shifts. These are also affected by the presence of the micro-sphere and 
thus, in what follows, we give a brief discussion of these quantities. 
Since the free space characteristics of coupled dipoles are well known, we 
concentrate on the effects produced by the micro-sphere. We normahze all 
frequencies to the free space decay rate denoted by ly^ with 

ro=:f^hn[n,.G<«>(r„r„«).Hj = ^. (17) 
2mc 3mc 

Thus quantities of interest will be 

(a) Micro-sphere dependent width y^K^ and shift /q^* ^f the dipole's 
frequency 
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(b) the radiative interaction 
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The last expressions of (18) and (19) follow from Eq.(ll) for radially oriented 

dipoles on z-axis and include the contributions from all possible modes. 

In what follows, we demonstrate giant enhancement in characteristics like 

due to the excitation of the whispering gallery modes. The modes associated with 
(16) are especially important. We begin by examining the extinction coefficient 
gexb which is a measure of the on-axis transmission for plane wave illumination. 
In terms of the coefficients A„ and B„ given by Eqs.(13), (14), gext is defined as 

Qe.,=--^V(2n + r,Rc(A„+B„). (20) 

The results of numerical computations for the a23i2,i67 mode are shown in Fig. 2. 
In order to assess the contribution of the WGM we have shown in Fig.2a the 
extinction coefficient Qext with (solid line) and without (dashed line) the resonant 
term. The contribution of the WGM aiaii.iev can easily be appreciated by noting 
the narrow (note the scale in the horizontal axis) resonance at (Oq . Moreover, the 
background has relatively insignificant contribution to gext- The results for 
and Q.J for equidistant dipoles (d^=d2=d) and for two normahzed distances, 
namely, d/X = 0.l and 1.0 are shown in Fig.2b. The effect of increasing distance 
of the dipoles from the surface is clear from a significant reduction of the 
enhancement factor K^^^ = iaax(\K^\) from 1349.81 to 154.47. Note that single 

dipole enhancement factor Kq^ has the same magnitude for the configuration 
under study. Analogous results for a^g^i mode are shown in Fig. 3. In contrast to 
the case of a23i2,i67 mode, now the background plays an important role. 
Moreover, the enhancement factors for the same distances, i.e., dl X= 0.1 and 
1.0, are now reduced to 3.4 and 0.07, respectively. The distortion of the curves in 
Fig 3b is due to a nearby WGM 344 j. We now turn to a direct comparison of the 
enhancement factors K^j as functions of the distance of the second dipole 

d2l X , while the first is kept on the surface {d^-0). The results are shown in 

Fig.4, which clearly demonsttates the advantages of using high Q modes of 
larger spheres. Note the resulting negative sign in Eq.(19) due to odd n for a39 1 
which is why the curves in Figs.3b and 4b are flipped upside down. We have also 
verified that very similar results are obtained if the dipole' s frequency is in the 

vicinity of another WGM, namely, ajj^ 157 . Our numerical results suggest a nice 
"single mode" approximation 

K^ + iQ.^ « KK^X^-iifO-mo)]-', ^^^^ 

+ iXlrf « iJcSTod [k - i{(0 - 0)a )]"^ , 
where K = C0qI1Q is the half width of the WGM and the sign + (-) is for the 
even (odd) mode. K^^ and K^^^ differ from each other by the sttengths of the 
mode functions for the WGM at and . 

We have thus shown that the radiative interaction between distant, otherwise 
noninteracting dipoles can be mediated via the excitation of whispering gallery 



modes of micro-spheres. For high Q modes and dipole frequency close to the 
mode frequency, our calculations suggest the apphcabihty of the approximation 
where the contributions from the surface dependent Green's function can be 
replaced by one arising only from a single WGM. Our results should be 
important in several contexts, for example, in the studies of energy fransfer 
between donor- acceptor molecule pairs. The transfer rate can be enhanced by 
several orders of magnitude, as it is proportional to the square of the dipole- 
dipole interaction, which, as shown, can reach giant values when mediated by 
whispering gallery modes of larger spheres. We further note that there exist 
several possibilities for the experimental realization of the theoretical predictions. 
One could adsorb the donor-acceptor pair on the surface of the micro-sphere. The 
other possibility, in the spirit of classic experiments of Ref.9, could be to let the 
micro-sphere fall between two parallel atomic beams. Finally we note that in the 
foregoing we have considered the case of weak radiative coupling between the 
atom and the vacuum of the electromagnetic field. The regime of strong coupling 
is interesting in its own right and leads to important consequences "'"'^ if such 
interactions are mediated by the WGM's. 

One of the authors (SDG) would like to thank the Department of Science and 
Technology, Government of India for supporting this work. 

APPENDIX A 

Quantum electtodynamics of the radiative interaction 

In this appendix we present a brief derivation of the radiative interaction. Let 
E(r,t) represent the quantized electtomagnetic field. In dipole approximation 
the interaction Hamiltonian of the atoms with the field is 
H,it) = -'£d„ ■Ei?„,t)\e„){g„ \e'''o'+H.C., (Al) 

a 

where we use the interaction picture and sum over the two atoms. The second 
order contribution to the wave function is 

I 'J' ) = j j "^^2^, ih )H, ih )| ) (A2) 

The initial state l*?*^*^^^ is |e^,g£,0), where |0) represents the vacuum state. The 
final state is | g ^ > 'O) hence the ttansition amphtude A will be 
1 r \ 

A = -—^dt,^ dt^{g^eBfi\H,{t,)H,{t^)\e^,gBfi). (A3) 

A calculation shows that the rate of change of the ttansition amphtude in the long 
time limit is given by 

dA 



-j^"\{d:.E{f^,r)d,.E{?„0))e-'''o'dr, 



° (A4) 



which on using the properties of the electromagnetic field vacuum can be written 
as 



The last line in (A5) follows from the definition of propagators in quantum 
electrodynamics in terms of the commutators of the electric field operators. Note 
the close relation of (A5) to Eq.(7), which we derived from classical arguments. 
The above argument clearly shows the relation of radiative interactions to the 
propagators. The argument leading to (A5) is quite general. Thus if the medium 
between the two atoms is something else, then the radiative interaction will be 
obtained by using the modified propagator. 
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Figure Captions 



Fig.l Schematic view of the micro-sphere of radius a with two similarly 
oriented dipoles on the z-axis, at distances of and from the surface of 

the micro-sphere. Refractive index of the sphere material is yfe = 1.46 , with 
that of outside medium ^Je^ = 1.0. 
Fig.2 (a) Extinction coefficient Q^^, and (b) dipole-dipole characteristics 
(solid line) and (dashed line) as functions of detuning cOq - (Oq for 
d/X = 0.l for a23i2,i67 mode. Both the curves in (b) are normahzed to the 
enhancement factor A'o^ = max(|A'^|) =1349.81. The corresponding curves 

for d/X = l.O wee identical except that now K^j =154.47. 

Fig.3 Same as in Fig.2 except for a^g i mode. The corresponding enhancement 
factors for distances d/ A,=0.l and 1.0 are now 3.4 and 0.07, respectively. 

Fig.4 Dipole-dipole enhancement factor K^^ in the vicinity of the cavity 
resonance as a function of t/j for (a). a23i2,i67 and for (b) aag i modes. The 
first dipole is on the surface. The dashed lines show the same when the 
resonant mode is suppressed. At d2/A. = 0, these values are 0.00394 and 
0.070, respectively. 




Fig. 1 G S Agarwal and S Dutta Gupta, 
"Giant radiative..." 
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Fig.2. G. S. Agarwal and S Dutta Gupta, "Giant radiative ..." 
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